Many bosons in a narrow resonance: the f?*-phase issue 
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The many-boson problem in presence of an asymptotically narrow Feshbach resonance is consid- 
ered. The energy spectrum is shown to be bounded from below in the limit of zero range interaction. 
This implies the promising possibility of achieving a strongly interacting bosonic phase in a dilute 
regime where the details of the actual interatomic forces are irrelevant. The integral relation between 
the energy and the one-body momentum distribution is derived. 
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The magnetic Feshbach Resonance (FR) technique al- 
lows the study of highly correlated and dilute many-body 
systems in which the two-body scattering cross section is 
arbitrarily large. This has led to ways of studying the 
cross-over between fermionic (BCS) and bosonic (BEC) 
superfluidity in the two spin-component Fer mi gas [l|| and 
to obtain the first evidence of Efimov states |2h12| as re- 
viewed in Ref. [13| . However, large atomic losses due to 
three-body recombinations prevent the observation of a 
thermalized resonant Bose gas. Nevertheless, the large 
scale separation between the scattering length (denoted 
by a) and the interaction range (denoted by b, of the 
order of few nanometers) suggests the existence of uni- 
versal low energy properties in this system flil [l5j . The 
prediction [ilj - Eol ] and the observation [2l|, l22j| of univer- 
sal four-body states are the first signatures of possible 
universal features for a many-boson system in presence 
of three-body Efimov states. The universal character of 
many-body properties at energies much smaller in abso- 
lute value than Eb = ti 2 / (Mb 2 ), should be revealed using 
a model where the interaction range is arbitrarily small. 
Unfortunately, structureless point-like interactions lead 
to the so-called Thomas collapse which makes the stabil- 
ity issue of the resonant Bose gas puzzling (23[. A well 
defined low energy model for the resonant Bose gas is 
thus needed. Eventually unlike the fermionic case, little 
is known about the resonant Bose gas which is now a 
fascinating issue at the cutting edge of the researches in 
ultracold atoms. 

Recently, it has been shown that three bosons in an 
asymptotically narrow FR experience an Efimov effect 
without Thomas collapse in the limit of zero-range inter- 
action (24[. In this system, the three-body parameter is 
only a function of the so-called width parameter R* and 
is independent from the details of the actual interatomic 
forces [24], |25[ . Moreover low energy states are predicted 
to be long lived 26] . Inspired by the three-body problem, 
asymptotically narrow FR are a relevant way to achieve 
a thermalized and dilute resonant many-boson system. 
In this letter it is shown that in such a regime the system 
has a well defined ground state, denoted as the i?*-phase, 
without Thomas collapse. The tail of the atomic mo- 
mentum distribution at high momentum k is shown 
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Furthermore, the energy E for bosons of mass M verifies 
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where (K mo1 ) is the mean translational kinetic energy of 
the molecular state and E tYap is the contribution coming 
from the interaction with a possible external potential. 
Equation |2) generalizes the Tan relation of Refs. 27[ to 
the case of an asymptotically narrow FR. 

The FR mechanism is based on the coherent coupling 
of atomic pairs (belonging to an open channel) with a 
molecular state (belonging to a closed channel). Us- 
ing an external magnetic field B permits the tuning of 
the molecular binding energy E mo \ with dE nlo \/dB = <5/i. 
As a result, the scattering length in the vicinity of the 
FR located at B = Bq can be arbitrarily fixed through 
the law a = Ob g [l- AB/{B-B )}. Here, AB is the 
width of the resonance and is the background scat- 
tering length. For a weak coherent coupling between 
the atomic pairs and the molecular state, the reso- 
nance is narrow (\5[iAB\ <C E b ) and the width param- 
eter R* = h 2 /(Mcib g 5[iAB) becomes essential in the de- 
scription of low energy properties. Asymptotically nar- 
row resonances where R* is very much lar ger than b, are 
characterized by the scattering amplitude [24| 



/ (fc) = -(l/o+ R*k 2 +ik)' 
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where the actual range of the interaction b is negligible. 

In what follows, the system is described by using a 
generic two-channel model which contains the FR mech- 
anism 28]. In this model, the molecular state (of mass 
2M) is structureless: a reasonable approximation for the 
description of mechanisms at energies much smaller than 
Eb. The operator a£ (or b\j creates an atomic (or a 
molecular) plane wave (r|k) = exp(ik • r) and obeys the 
usual commutation rules for bosons. The kinetic energy 
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of a particle of wave vector k is denoted by = , so 
that in absence of an external potential the free Hamil- 
tonian for atoms and molecules is 



Hn = 
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The interchannel coupling converts a molecule into a 
pair of atoms via the operator V and vice-versa via the 
operator V' . It is modeled by a constant amplitude 
A (chosen real positive) and a generic cut-off function 

Xe (fc)=exp(-^) (withe>0) ; 
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The Hamiltonian is thus Hq + V + and the scattering 
amplitude denoted by fe(k) at energy E — E^ is given by 
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where the variable z encapsulates the usual prescription 
z = E + i0 + for E > and z = E otherwise. For vanish- 
ing values of e, identification of / c (fco) with Eq. ([3]) gives 
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In what follows, only states with asymptotically small 
values of e are considered while a and R* are kept fixed. 
In this zero range limit, Eq. ([7]) shows that E mo \ has 
an arbitrarily large positive value. Thus the molecular 
state is occupied only through virtual processes. As a 
consequence of the interchannel coupling, an eigenstate 
for a system composed of N atoms is in general a coherent 
superposition of states containing to [N/2] molecules 



[N/2] 
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where [ip™) is a state with m molecules and N — 2m 
atoms. Projection of the stationary Schrodinger's equa- 
tion at energy E in the subspace with to molecules 
(1 < to < [N/2]) yields the recursive equation 
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The recursion begins for to — by the equation 
|^) = l*i„c)+Go(z)VlV>!>, 



(9) 
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where Gq(z) = (z — i?o) _1 is the resolvent and l^mc) is 
the possible incoming wave for a scattering process when 
E > whereas |<l>inc) = when E < 0. The recursion for 
to = 1 plays a central role in the rest of this letter. From 
Eqs. © and {JTOj) it takes the form 

[E-Ho- V^Go(z)VM) = V^|$ inc ) + V\^). (11) 



In the subsequent lines, the momentum representation is 
used and the hyper wavevector (k) p = (ki, k 2 . . . , k p ) [or 
(K) m = (Ki, K 2 . . . , K m )] denotes the wavevectors of p 
particles (or of to molecules) . The m-molecule state can 
then be written as 



We 
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In the following, a superscript i in (k) p means that the 
wavevector k^ is withdrawn from (k) p . For instance 

= (k 3 ,k 4 , . . . kjv). The ket \(k) p ] (K) m ) describes 
the state where the atom i (1 < i < p) is in the plane 
wave of wave vector k^ and the molecule j (1 < j < m) 
is in the plane wave of momentum Kj. Furthermore 

kjj = *o 3 (or Kij = ki + kj ) denotes the relative (or 
the total) wave vector for the pair 

As in the two-body states, for an arbitrarily large rel- 
ative wavevector of the pair (ij), it is expected that 
in the presence of to + 1 molecules (i.e. \ip™ +1 ) ^ 0) 
((k) v ; (K) m \il;™) has the asymptotic behavior Xe(hj)/k^j 
29]. This behavior is obtained by considering Eq. (j9j 
divided by E mo \ = 0(1/ e) in the zero range limit; the 
only non vanishing contribution in the right hand side 
of the resulting equation is ((k) p ; (K) m \V^\ip™~ 1 ). This 
enables the TO-molecule wave function to be expressed in 
terms of the (to — l)-molecule wave function: 
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In the center of mass frame, can be factorized by 
a ket of N — 2 particles, which is denoted by \D) in the 
zero range limit: 

lim ((*)$; Ki 2 |#> n 

£ ^ 0+ - =8(y k i +K 12 )((fc)] v 2 |J). (14) 
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Similarly, the possible incoming state is factorized as 

: ) . The equation 



((k) N \^nc) ^ (27r) 3 S(j:Zi^)((k)N\^ 
for \D) is deduced from Eqs. (f5]> . ffTTjl and ([13]) and the 

identity lim^oo fc?-((fe)^|£>) = lim^-^ kf 2 ((k)%\D) 
which comes from the exchange symmetry of the two 
molecules described by \ip 2 ) for {i, j} ^ {1,2} [301 ] : 
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where fc co1 = V ME col /h is the collisional wave num- 
ber for the pair (1,2) defined from the energy 
E co1 = E — £^ =3 £ k „ - E Kl2 /2. Surprisingly Eq. © 
supports solutions where lim e _j.o IV^™) = f° r m > 2 and 
thus limj^oo l${{k)%\D) = for {i,j} $ {1,2} 

Without loss of generality, the energy spectrum of 
Eq. (IT51) is shown to be bounded from below in the strict 
resonant limit \a\ = oo. To show this, we use the fol- 
lowing inequality verified by negative energy solutions of 
Eq. (HSJ: 
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where q = \f—ME/fi is the binding wavenumber. Then, 
we consider a hypothetical solution of Eq. (fT5|) having an 
arbitrarily large and negative energy (q —¥ oo). Assuming 

that ((fc)^J|D) is a bounded function, integration in the 
domain where ki2 is smaller or of the order of q in the 
right hand side of Eq. (p~6|) gives a finite result in the 
limit where q is arbitrarily large. In the domain where 
ki2, is much larger than q, the high momentum behavior 

of ((fc)^|_D) [32j | gives rise to a vanishing contribution 
after integration in Eq. (|16j) . As for any fixed value of 

(&)]y the left hand side of Eq. ([T6]) is bounded from above 
by a finite quantity, gi?* cannot be arbitrarily large. The 
desired result is thus obtained by contradiction. 

From now the Effective Range Model (ERM) is used to 
derive Eqs. (JIJ and ©. To avoid any confusion with the 
two-channel model, the eigenstates in the ERM are un- 
derlined (as for instance). The ERM is a one-channel 
and zero-range approach which supports the scattering 
amplitude in Eq. © whereas the molecular field is hid- 
den. The main idea of the ERM is to filter asymptotically 
(i.e at interparticle distances much larger than b) correct 
wavefunctions among the set of singular wavefunctions 
satisfying a Schrodinger's equation for point-like interac- 
tions. This filtering process can be performed by using 
a pseudo-potential. Following Refs. [33j, |34|, the pseudo- 
potential of the ERM is defined by introducing a Gaus- 
sian shape (r\8 e ) tending to the Dirac distribution in the 
limit where the short range parameter e goes to zero. In 
the momentum representation (k|<5 £ ) = Xe(k) i.e. it coin- 
cides with the interchannel coupling function of the two- 
channel model. For a many-body system, a short range 
parameter and a pseudo-potential V^ 3 is introduced 
for each interacting pair of particles {i,j): 
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where the ket \S ei ) is associated with the function 6 ei . 
for the relative coordinates of the pair. Instead of intro- 
ducing the pseudo-potential, the filtering process of the 



ERM can be defined through the notion of contact con- 
dition. To see this, one can notice that the action of the 
pseudo-potential V^ lJ on an eigenstate |$} gives rise to a 
(5-source term for each interacting pair in the stationary 
Schrodinger's equation ; 
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Here, |^4^' 7 ' ) ) depends on the set of short range pa- 
rameters {e pq } such that (p,q) ^ {i,j) and verifies 

\A%) = - V^721 im ^o+ M<M*» SI- The ERM 
can then be defined from Eq. ([T8"l) together with a bound- 
ary condition at the contact of each interacting pair {i, j) 



lim B U] {8 U] m =0, 
where B t is the boundary operator defined by 
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From Eq. (fT8")) one finds that two non degenerate eigen- 
states of the ERM, say |\S[) (of energy Ey) and |$) (of 
energy E$ ^ Ey), are not mutually orthogonal: 
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Ant? 
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E,i, — Ea 
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Interestingly, a modified scalar product denoted by (-|-)o 
which makes the ERM self-adjoint can be constructed by 
using Eqs. p^|) . (j2H|) and (|2"T |) : 



(<H*)o = nm 

{e}->0 



i<3 



(22) 



where lim^^g denotes the limit in which all the short 
range parameters tend to zero. Equation (|22[) generalizes 
the modified scalar product introduced in Ref. (36l | for 
two interacting particles. The equivalence of this modi- 
fied scalar product with the usual scalar product of the 
two-channel model has been already depicted at the two- 
body level in Ref. [28[ . This mapping can be generalized 
in the A-body case as follows. In the zero-range limit, 
analogously to Eq. (fT3"|). the hidden m- molecule wave 
function is revealed in the state |\&) by the existence of a 
k~j 2 behavior in the large momentum limit for m distinct 
pairs {ij). To be more explicit, one introduces the state 
|^ttj),...(Jn),(pg)^ j^j corresponc iing to the contact of the 
two particles in the distinct pairs {ij),. ■ ■ {In), {pq). It is 
defined iteratively from the state \A^) by 



({k)%- lA ;K lJ ...K ln \A { * 



(ij),...{ln)\ 



Xt pq {kpq) 



k' 2 



x {{kyt tm ^i---^ln\^ ), -" iln),< ^ ) )- (23) 
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In the scalar product the contact for m distinct 

pairs (ij), ■ ■ ■ (In) manifests itself by the monom 



(-1)' 



16tt 2 



( 27r )3m/2 V 



..e ln (A^- {ln) \A^- {ln) ). (24) 



The number of monoms of the same type as in Eq. ([24]) 
is equal to = N\/[2 m m\(N - 2m)!] i.e. the number 
of collections of m distinct pairs taken from the set of 
N particles. Applying all of the boundary operators to 
the right hand side of Eq. (|22|) . and identifying the result 
with the scalar product in the two-channel model, one 
finds the following mapping in the zero-range limit: 



lim \A 
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and lim{ e }^ 1 .3D = lim e _>o+ IV'e)- Equation ([15]) is recov- 
ered using this mapping and applying the contact con- 
dition in Eq. (TH)1) on the state For an eigenstate 
normalized with Eq. ([22]) i.e (^_\3i)o — lj the one-particle 
momentum distribution (taking into account the contri- 
butions of the m-molecule states) is rt k = C3{.l a k a k|.3!.)o- 
This definition yields the correct normalization i.e. 



d 3 k 



r n k = N — 2N*> 



N" 



[N/2] 
m— 1 



_ (26) 

where iV mo1 is the mean number of molecules [38| . The 
ERM gives c 4 = 87rY mol /i?* in Eq. (U): a result which 
also can be deduced from the Hellmann-Feynman the- 
orem [39|, l4pj . The expression of the coefficient c§ in 
Eq. ([T]) depends on the possible interaction of the parti- 
cles with a trapping potential V tIllp (r): 
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Equation ([2]) is obtained using the identity 
C3il^* J l.3Do = 0, which is valid provided that the 
state \^_) verifies the contact condition in Eq. (|19|) for all 



pairs Thus, applying lim {£} ^ (l\. l<:j B tt ^j on the 

mean value of the Hamiltonian of the ERM evaluated 
for the normalized eigenstate yields 



E - £ trap = t£| 



d 3 k h 2 k 2 

(27t) 3 "2AT 



a k ak|*)o, 



(28) 



where £: tra P = (*| ^„ = i ^ trap (r„)|*) includes the con- 
tribution of the molecular state. The regularization of the 
kinetic energy in the right hand side of Eq. (f28| through 
the modified scalar product permits to recover Eq. @ 
with (X mo1 ) = (*| f 6 k S k |*)/2. 

To conclude, the main result of this letter is the sta- 
bility of the resonant dilute many-boson system in an 
asymptotically narrow FR. The system is described in 
terms of a two-channel model of the FR. The mapping 
with the ERM permits to obtain the integral relation 
satisfied by the energy and the one-particle momentum 
distribution. In the ERM, the hidden molecular field is 
described through pairwise singularities in the iV-particlc 
wavefunction. The absence of Thomas collapse for three 
bosons in an asymptotically narrow FR is at the heart 
of the stability issue of this many-boson problem. The 
attractive effective force between one boson and one pair 
of bosons close together, is screened for interparticle dis- 
tances smaller than the width parameter R* [2a |. Be- 
cause of this feature one expects the system to be orga- 
nized in pairs of bosons with a large population of the 
molecular state and a mean intermolecular distance of 
the order of R*. The liquid or crystalline nature of the 
i?*-phase as a function of the particle density is an open 
fascinating issue. It is worth pointing out that Eq. ([3]) 
is valid in a small interval of magnetic field such that 
\B — Bq\ <C |AB|. Otherwise, off-resonant effects must be 
taken into account 4l| . Ultra narrow FR (where R* 3> b) 
have been already found in many atomic species [42| and 
new techniques allowing a fine tuning of the magnetic 
field are promising for the observation of the i?*-phase 
0. 

Stimulating discussions with Yvan Castin, Tom Kris- 
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